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The dissipative conductance of an array of mesoscopic rings, subject to an a.c. 
magnetic flux is investigated. The magneto-conductance may change sign between 
canonical and grand-canonical statistical ensembles, as function of the inelastic level 
broadening and as function of the temperature. Differences between canonical and 
grand-canonical ensembles persist up to temperature of the order of the Thouless 
energy. 
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Many of the idiosyncrasies of mesoscopic systems may be found in the response of small 
metallic rings to an a.c. bias. Previous works in this field have stressed the role of intralevel 
(diagonal) and interlevel (off-diagonal) transitions as well as the fact that the overall 
physics may depend on the type of external bias [^] (coupling to external reservoirs, ex- 
ternally applied electromagnetic fields, etc.). Several approaches have been taken to derive 
quantitatively the average multichannel conductance of disordered mesoscopic conductors 
(in the diffusive regime). We mention here two. One is based on impurity Green's func- 
tion techniques. In the context of Aharonov-Bohm cylinders, the leading order quantum 
corrections to the Drude conductance has been first calculated in Ref. Among other 
interesting results of this theory, it also predicts small positive magneto-conductance (MC), 
in conjunction with other weak localization studies and a semiclassical picture. Another 
approach that has been extensively explored is based on the observation that in disordered 
conductors, spectral properties of single electron levels (within a not-too-large energy range) 
are satisfactorily described by random matrix theories PJ^l- Remarkably enough, results of 
such a random matrix approach for the a.c. absorption (hence the real part of the conduc- 
tance, g), differ sharply from those of the former approach. For one thing, a negative weak 
field MC is predicted. To the best of our knowledge these differences hitherto have passed 
practically unnoticed. 

The purpose of our systematic study, reported here, is two fold. By extending previous 
analyses we have attempted at bridging the apparent differences among the various ap- 
proaches alluded to above, presenting a coherent picture of the quantum a.c. conductance 
and its dependence on important parameters of the system. At the same time we have inves- 
tigated differences between the canonical and the grand-canonical statistical ensembles, (CE) 
and (GCE) respectively, which have proven to play a crucial role in the thermodynamics of 
such systems 0. 

The central point of our discussion consists of the observation that two major factors 
conspire to determine the behavior of our system (where, throughout most of the present 
analysis the sign of the MC is taken to be a signature of this behaviour). The first is the 
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dimensionless parameter A/7, where A is the average level spacing at the Fermi energy and 
7 is the relaxation rate (see below) which determines the coupling of the electronic system 
to the eternal degrees of freedom, hence the broadening of the individual single electrons 
levels. Technically, for A/7 < 1, (> 1) one employs perturbative i.e. Green's function, (non 
perturbative, e.g random matrix theory) techniques. But also the physics differs sharply 
between the mesoscopic continuous spectrum limit and the microscopic discrete spectrum 
limit. The second important factor is the type of the statistical ensemble employed. We 
argue that even when dynamical response is concerned different terms may contribute to 
the conductance depending on whether GCE or CE is employed. Consequently in the 
non-perturbative limit the MC may differ by sign and order of magnitude between the 
two ensembles. We also identify such differences, albeit small (of order (A/7)^), in the 
perturbative regime. 

Our main results are summarized at the end of this paper. For the sake of specificity we 
consider here quasi one-dimensional cylinders [L^ > I, where is the perimeter and / - 
the elastic mean free path. The inelastic mean free path, l^, is larger than all linear scales 
of the cylinder). Our analysis is extendible to other geometries (including quantum dots) 
where by and large, we anticipate our qualitative predictions to hold. 

Within linear response theory, we consider a small amplitude flux component at frequency 
uj superimposed on a static Aharonov-Bohm flux, $. The average conductance of a cylinder, 
g, is related to the conductivity by a factor S/L^, S being the cross section area for the 
current density @|. Our starting point follows that of Ref. [^. The density matrix p satisfles 
the kinetic equation 

^ + ^[H,p] = -l{p-pQE). (1) 
Here H is the Hamiltonian; 7 is the relaxation rate towards an instantaneous equilibrium 



state described by pqe A straightforward apphcation of linear response, (cf. Ref. |^), 
results in the following expression for the conductance g 
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where e„ = en{^) and |r;, >= |?t,($) > denote the exact single electron eigenenergies and 
eigenstates; /„ = — is the Fermi-Dirac function [|ri|], and J = —dH/d^ 

is the current operator. We shall denote the respective impurity averaged terms on r.h.s. 
of Eq. by gi,gn,giii and gjy- We note that gj is basically the flux derivative of the 
persistent current, yielding a purely imaginary contribution to g, which has been evaluated 
both within the GCE and CE [O. As a preliminary we note that for the canonical ensemble 



gii cancels against gm at zero temperature. This implies that in a canonical situation there 
is no diagonal relaxation at zero temperature (the occupation probability of any given level 
is flux independent). Hereafter we consider the frequency range ut <^ 1, where r is the 
elastic mean free time. 

Diagrammatic approach (valid for 7 ^ A). Within the grand canonical ensemble 
/x($) = 71 is a constant, hence gu vanishes. The remaining contributions of gni and giv 
yield |] 

" 1 A /o/ — 7V, I 
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where Fi(x,$) = ia;sinhx/(coshx— cos47r$/$o) and (70 = ^"^"^ is the Drude conductance. 

Here = hD/L^ is the Thouless energy, D is the diffusion constant. The second term in 

Eq. (^, first derived in Ref. |^ is the A/7 weak localization correction to the Drude result 

|T^. We recall that indeed systematic diagrammatic theory analysis may be presented as 



an expansion in the small parameter A/7. The temperature, dependence of g^'"^ comes 
mainly from the dependence of 7 on T. The real part of the conductance, ^g'-^^^, gives rise 
to positive MC at small flux (see inset to Fig. 1). 

Within the canonical ensemble the number of electrons in each ring remains unchanged 
as the flux is varied, which implies $ dependence of /x. To evaluate the term g^j we neglect 



the $ dependence of / (giving rise to higher order terms in A/7). We obtain at T <^ 7 |14 



^f/ = ^o^^-F2(#,$), (4) 



where F2{x, $) = x'^^ ^^^^'^'' J^^^^''^\ The temperature dependence of g'"^ is described below. 
We note the correspondence between gff" and the typical single level current: = 



-47r{AEc/<l>l)F2{^A/{'KEc), $) |12|,|T5l. To evaluate the last two terms of Eq. (g), gfjj+gf^ 
we expand fn about the flux and impurity average chemical potential JI. The result can be 
written as [|T^: 

gffi + 9?v^ = + ^\viuJ,^-r^{0M, (5) 



—lUJ 



with 



Evaluation of rj necessitated the calculation of a considerable number of diagrams with 
numerous energy integration ranges. A straightforward but lengthy calculations give 



_ 2 
ri{uj,(l>) = go — 



[$ ^ 0]. 
(6) 



Hereafter we shall discuss the dissipative part of g. We note that the contributions Eqs. 
(I) and to the MC are comparable and of the same sign. Indeed $) are oscillatory 
functions of $ with a period $o/2 and an amplitude of order unity (for x < 1). We stress that 
all terms except g'^^^ (i.e. all terms particular to the CE), contain an explicit temperature 
dependence. We find that their contribution is roughly a constant ~ (70(^/7)^ for T <^ 'j, 
then decays as a power law ~ goA^ / {iT) ioi •y T <^ Ec and finally vanishes exponentially 



~ goAym expi-^T/E,) for E, « T 0. 

The flux dependence of the canonical terms ^{g'"^ — g'-^'^^'j is depicted in Fig. |l] for few 
values of T and 7. We note that the canonical terms exhibit flux dependence opposite to 
the commonly accepted GCE weak localization behavior (thus showing the tendency to the 
negative MC). This dependence is eventually suppressed at temperatures T > Ec. 

Non— perturbative regime. For 7 < A one is not able to employ perturbative tech- 



niques. We thus adopt another approach following the treatment of Refs. [§[ and [|T8| of 
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the a.c. absorption in an applied electric field. Certain important modifications are due. 

Within the GCE the real part of the diagonal contribution gffj^ {gff^ = 0) is given by 

GCl _ 



-,GCE _ 2e^ <ilm>H rj^j^ j^^^ equality may be written as (cf. Eq. (|)), 
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To obtain Eq.(|3) we have employed some of the results of Ref. 1 12, 15, 19 1- This evidently 



leads to a positive MC, with an amplitude 2gQ (^^^ t^i;^ I|20|. For <^ 7 we recover the I/7 
divergence discussed previously for intralevel absorption [|l],H|. Unlike in the perturbative 
regime, here this is the dominant contribution. 

To calculate the off-diagonal contributions we first note that we may perform separately 
the ensemble average of | < n\J\m > p and the terms that depend on the eigenvalues e„, 
(see Eq.(0)). The former being averaged over an energy interval Ec is replaced by a constant, 
whose value is determined to be goA'^/n, compatible with the requirement that g ^ go for 
7 > A PJ2T|J21]. We thus may write giv as 



„GCE _ 90 
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where (. . .)gc£; refers to averaging under grand canonical conditions (e.g. assuming /i to 
be uniformly distributed over an energy interval much larger than A). The quantity in the 
brackets is i?(e, $) - the level pair correlation function. It may be shown that , to leading 
order in x = vre/A |2|j23 
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with V = (^TT^'^I'^ ^Ef-I A ($/$o)- The leading contribution to g^y^ are respectively (x to be 
replaced by z = ixioj + i7)/A [|^): gozijT + 2i(ln2; + F — 2))/6 in orthogonal; g^zi^iX^nv + 
y/TTz/v)/Q in intermediate; and goz{z — 2i)/3 in unitary cases, where F is Euler's constant. 
This off-diagonal term is marked by a negative MC of amplitude go{'y/A) ln(7/A) (at u = 0). 



We note, though, that this off-diagonal term is only a small correction on top of the diagonal 
contribution, (the latter being of order of A/7 when = 0, cf Eq. (|^). The total MC 
(including the dominant A/7 diagonal contribution Eq. (|^) gives rise to a positive MC. 
There is practically no explicit T dependence in the GCE case. 

The situation is markedly different for the CE. As was noted above, the two diagonal 
contributions to g*-"'^ are offset at T = {g^^^g = gff + g^fj = 0, for T = 0). This offset 
is highly temperature dependent. The difference Ig^j^g — gdSi.g\ decreases significantly at 
T ^ A, and continues to decrease further as until at T > Ec differences between two 
ensembles disappear. Thus the diagonal contribution, giving rise to positive MC, dominates 
at sufficiently high T see Fig. ^. 

As for the off-diagonal contribution, the average in Eq. (||) is to be replaced by (. . ■)ce, 
implying that particle number, rather than chemical potential fi, is a uniformly distributed 
random parameter. We note that (. . .)ce ~ (• • ■)gce Making further the approxi- 
mation {n — m)A ^ Asigne + e, we obtain (cf. Eq. (|)) 

= + «i /" i.^^^ (10) 

in J-oo e — — «7 |e| 

Note the extra |e|~^ factor in Eq.(|TUp. Consequently, the sensitivity to small magnetic 
flux (arising due to change in a level statistic), is more pronounced in the CE than in the 
GCE. We obtain for the r.h.s. of Eq. ([TO|) g^n'^ /Q in orthogonal; gQ^/^T{^TZ + 2i In z) / (Qu) in 
intermediate; and go{nz + 2zln2;)/3 in unitary regimes. At T = 0,ct; = this yields for 
the real part of conductance, goir'^/G for orthogonal and (7o(27r/3)(7/A) ln(A/7) for unitary 
cases. At T = 0,7 = we obtain g^n'^/Q and gQin"^ /?>)uj respectively [|^. This leads to 
a large negative MC of amplitude g^ (to be compared with the small corresponding off- 
diagonal contribution in the GCE case). Numerical results for the flux dependence of g^^ 
are shown in Fig.^. Our numerical data are compatible with our analytical results, showing 
for the CE a transition from negative to positive MC as T increased. More details of our 
study will be published elsewhere. 

In summary, we have reported here a comprehensive study of the conductance of meso- 

7 



scopic rings in various regimes of the relevant parameters. In particular we have noted that: 
(1) Within the CE the MC may change sign as function of the parameter A/7, going from 
the discrete level to the continuum limit. (2) The incipient transition may be seen from a 
perturbation expansion in powers of (A/7). The (A/7)^ term possesses a sign opposite to 
the (A/7) term and it eventually takes over as (A/7) approaches unity. (3) The diagonal 
term with the CE is highly temperature sensitive. In the A/7 > 1 regime this gives rise to 
a sign change of the MC as a function of temperature. Within the GCE the corresponding 
term is by far more robust to temperature, rendering the MC always positive. (4) Differ- 
ences between the two ensembles, accentuated at low temperatures, diminish as a power law 
with temperature up to T Ec- (5) We have found contributions to the conductance which 
mainly depend on elastic scattering (e.g Eq. (p!OD), while others are dramatically dependent 
on the inelastic rate (e.g. the diagonal contribution Eq. (0)). This sheds light on previous 
arguments concerning the relative importance of elastic and inelastic processes in determin- 
ing g [|1],|2[ . We finally note that within present technology, one is not too far from achieving 
the limit 7 < A in experiment P? . 
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FIGURES 

FIG. 1. Flux dependence of the canonical terms ^{g^^ — g'^^^) (in units of /%) for 
EJA = 2Q,u = 0. (Solid line) T = 0,7 = 5A; (dashed line) T = 0,7 = lOA; (doted line) 
T = 10A,7 = 5A. Inset: the GCE conductance (solid line) 7 = 5A; (dashed line) 7 = lOA. 

FIG. 2. Flux dependence of the canonical average conductance (in units of e^/fi), obtained 
from numerical simulations on the Anderson model for a 3D ring of dimensions 30 x 4 x 4. Here 
E'c/A = 7/27r,a; = 0, 7 = .3A. The numerical results are averaged over energy range and over two 
disorder configurations. The data obtained for the CE at T = Ec are identical with the GCE. 
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